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ABSTRACT 

We present a study of protoclusters associated with high redshift radio galaxies. We imaged MRC 1017-220 (z = 1.77) and MRC 0156- 
252 (z = 2.02) using the near-infrared wide-field (7.5' x7. 5') imager VLT/HAWK-I in the Y, H and Ks bands. We present the first deep 
K-band galaxy number counts within a large area (~200 arcmin 2 ). We then develop a purely near-infrared colour selection technique 
to isolate galaxies at 1.6 < z < 3 that may be associated with the two targets, dividing them into (i) red passively evolving or dusty 
star-forming galaxies or (ii) blue/star-formation dominated galaxies with little or no dust. Both targeted fields show an excess of star- 
forming galaxies with respect to control fields. No clear overdensity of red galaxies is detected in the surroundings of MRC 1017-220 
although the spatial distribution of the red galaxies resembles a filament-like structure within which the radio galaxy is embedded. In 
contrast, a significant overdensity of red galaxies is detected in the field of MRC 0156-252, ranging from a factor of ~ 2 - 3 times the 
field density at large scales (2.5 Mpc, angular distance) up to a factor of ~ 3 - 4 times the field density within a 1 Mpc radius of the 
radio galaxy. Half of these red galaxies have colours consistent with red sequence models at z ~ 2, with a large fraction being bright 
(Ks < 21.5, i.e. massive). In addition, we also find a small group of galaxies within 5" of MRC 0156-252 suggesting that the radio 
galaxy has multiple companions within ~ 50 kpc. We conclude that the field of MRC0156-252 shows many remarkable similarities 
with the well-studied protocluster surrounding PKS1 138-262 (z=2.16) suggesting that MRC 0156-252 is associated with a galaxy 
protocluster at z ~ 2. 

Key words, large scale structure - galaxies: clusters: general - galaxies: evolution - galaxies: high redshift - galaxies: individuals 
(MRC 1017-220; MRC 0156-252) 
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'■O 1. Introduction sive and star-forming galaxies occurred in clusters, one needs to 

^— ' ■ . , , • , • , study overdensities at higher redshifts. 
t-H . It is well established that the evolution or galaxies strongly 

> depends on environment. In the nearby Universe, the high- 

est density regions (e.g., the cores of local galaxy clusters) However, the number of known galaxy clusters at high red- 

^ • are dominated by red, early-type galaxies with the fraction of shlfts 6 > 18 smalL Indeed ' marching for higher redshift clus- 

^ ' blue galaxies with on-going star-formation significantly smaller ters ra P ldlv becomes challenging using classical detection meth- 

«J ! than lower den sity regions (Dressier 1980; Tanaka et al. 2005; ods such as usln 8 red-sequence algorithms aimed at identifying 

iPostmanetal . 2005; Balogh et al. 2007; Poggianti et al. 200<C overdensities ^ladders & Yee| | 2000t | Andreon et al. | 

The spectra of these red passively evolving galaxies show a R^g ^J^^ ^ ^^.S ^g^- 
characteristic break at 4000A — i.e., the light from old stel 



cluster medium dStanford et al.ll2006t iRosati et al.ll200l l2009h . 
Until recently, the two highest redshift clusters were discov- 



lar populations is more prominent than that emitted by younger , , ,/ v vhuv™ munnio * , 

, * *L, i i- . . . - i ered through X-rays: XMMXCS J2215. 9-1738 at z = 1.457 

stars. These galaxies he on a tight red sequence in colour- , . , ,„ ° • „ ^ ■ , • ■ ■ ■ , 

... „. ,. , , . . , (with 17 spectroscopically confirmed members within the clus- 

magmtude diagrams. Studies have shown that this red sequence . . , r ,. \ „y . ' ^ . i \ ~„„a 1TT7L i j 

. , , „ i j - i r i ■ . . . . ter vinal radius; Stanford et al. 2006; Hilton et al. 2007) and 

is already well populated in clusters of galaxies out to high red- v . T ^ . : ~ n , .' . .. 

u-ft / 1 c c, e a t i ™c oaa2 iv / r • * i h nn/onnn XMMU J2235. 3-2557 a t z - 1.39 (w ith 34 spectroscopically 

shifts (z ~ 1.5; Stanford et al. 2005, 2006; Mei et al. 2006, 2009; ^ , , „ . . „ . . . . r J r 

rr-rr — ilUnnoLli^ i t llUnnnk t a t a i. n.- confirmed members; Rosati et al. 2009). Galaxies in the core of 

Lidman et al. 2008; Kurk et al. 2009). To understand when this ^ . , . ' — z — ~. rr, „ , , , , 

— ; — , . f — ~r , this second cluster already he on a well defined and tight red se- 

red sequence appeared, i.e., when the segregation between pas- rr-r: — ^ . l UrvXok t,, v , . , , ° 

2 quence (Lidman et al. 2008). These X-ray selected clusters were 

* Based on observations obtained at the European Southern however recently supersceded by the discovery of a galaxy clus- 

Observatory using the Very Large Telescope on Cerro Paranal through ter at z = 1.62, C1G J0218-0510, using photometric redshifts. 

ESO programs 081A-0673(A) and 083. A-023 1(A). " It is located in the Subaru/XMM-Newton deep field and has 15 
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confir med members to date dTanaka et al]l2010t iPapovich et al.l 
l2010h . 

One of the most efficient methods to search for galaxy clus- 
ters at even higher redshifts (z > 1 .5) is to look in the vicinity 
of high-redshift radio galaxies (HzRGs hereafter). These galax- 
ies are a mong the most massive galaxie s in the Universe (M> 
10 n M o ; iRocca-Volmerange et al] 12004 ISevmouret all [2007 ) 
and are g ood tracers of high densit y regions in the early uni- 
verse (see Mil ev & De Br euck 2008, for a review on HzRGs and 
their surroundings). Narrow line emitter surveys of HzRGs at 
2 < z < 5 show that they are often located in overdense regions, 
designated as 'protoclusters' that are likely to b e the progenitors 
of the present day massive groups and clusters ( Venemans et alJ 
l2002Ll2005Ll2007HPentericci et al.ll2000HKurk et alJl2004 . The 
narrow line emitters are, however, only a low-mass subset of 
the general population of UV-selected , star-forming galaxies 



ral p( 

dMilev et al.l l2Q04t lOverzier et~al1l2008l) and likely do not rep- 
resent the majority of the total stellar mass. 

Studies have tried to identify more massive cluster galax- 
ies associated with these radio-loud sources by looking for the 
red passively evolving ga laxies that may be populating the cores 
of high redshift clusters. Bes t et al.l d2003l) observed, for exam- 
ple, the environment of powerful radio-loud sources at z ~ 1 .6 
and found overdensities of red galaxies (R - K > 4) on two 
scales around the AGN: a pronounced central peak (within 
150 kpc) and weaker excesses betwe en 1 and 1.5 Mpc radius. 
More recently, iGalametz et alJ ((2009) studied the environment 
of 7C 1756+6520, a radio galaxy at z = 1.416, and found an 
excess of red sources (passive, early-type galaxy candidates at 
Z > 1 .4) within 2 Mpc of the HzRG A galaxy cluster associated 
with 7C 1756+6520 h as since been spectroscopically confirmed 
(IGalametz et alJl2010h . At higher redshifts (z > 2), studies have 
searched for red evolved galaxies by bracketing the redshifted 
4000A break with near-infrared filters. Kaiisawa et al. (2006) 
explored the environments of six HzRGs at z ~ 2.5 and isolated 
the evolved galaxy p opulation at z > 2 usin g purely near-infrared 
(JHKs) colour cuts. Kodam a~et al.l J2007) used the same near- 
infrared criteria to select protocluster mem ber candidates i n the 
field of HzRGs at 2 < z < 3 (see also IZirm etal]|2008l for 
a near-infrared study of a forming red sequence in a protoclus- 
ter at z — 2.16). They both found that some of their targeted 
fields contained overdensities by a factor of 2 - 3 compared to 
blank fields. Recent spectroscopic follow-up of these red sources 
has been conducted in a couple of these protocluster fields by 
iDohertv et alJ d2010l) . They confirm two red galaxies associated 
with PKS 1138-262 at z = 2.16, a dusty star-forming galaxy 
and an evolved galaxy with little ongoing star formation. These 
HzRG companions h ave an estimated mass of 4 - 6 x 1O U M . 
IDohertv e t al. (2010j) also confirm that a pure near-infrared cri- 
terion is efficient at selecting high redshift galaxies — e.g. 56% 
of their ///^-selected galaxies with spectroscopic redshift fall 
at 2.3 < z < 3.1. However, the low success rate of their spec- 
troscopic campaign confirms however the challenge in deriving 
redshifts for passively evolving galaxies whose spectra do not 
show prominent and easily identifiable emission lines. 

Building on these previous studies of individual HzRG 
fields, we defined a first uniformly selected sample of the most 
powerful radio galaxies in the pivotal redshift range 1 .7 < z < 
2.6 where we expect the cluster galaxies to start settling on the 
red sequence. Our method is to select potential cluster mem- 
bers (in particular, evolved passive galaxies) using near-infrared 
colour cuts. We observed our sample with the High Acuity Wide 
field K-band Imager (HAWK-I ; iPirard et alJl2004t ICasali et alJ 
2006; iKissler-Patig et alJ l2008h on the Very Large Telescope 



(VLT) in 2008 and 2009 in a set of three filters (YHKs or JHKs 
depending on the redshift of the targeted HzRG). This paper re- 
ports the results on our two lowest redshift targets. A companion 
paper presents results on the higher redshift targets (Hatch et 
al., 2010). 

We design a news near-infrared criterion to isolate galax- 
ies at 1.6 < z < 3. We apply this selection technique to 
study the galaxy population in the vicinity of the lowest redshift 
HzRGs of the HAWK-I sample: MRC 1017-220 (z = 1.77) and 
MRC 0156-252 (z = 2.02). We first describe in §2 the multi- 
wavelength data available for the targets and four control fields 
(including two sub-fields of GOODS South). The extraction of 
the source catalogues and galaxy number counts are described 
in §3 and §4. Section 5 presents the colour-colour selection 
technique we developed to isolate cluster member candidates at 
Z > 1.6. The properties of these cluster candidates, such as over- 
densities and spatial distribution are detailed in §6. Section 7 
summarizes our results. 

We assume a ACDM cosmology with Hq = 70 km s _1 
Mpc -1 , Q m = 0.3 and Q.\ = 0.7. All magnitudes are expressed 
in the AB photometric system unless stated otherwise. 



2. The data 

2.1. The targets 

We present a study of the two lowest redshifts targets of the 
HAWK-I sample: MRC 1017-220 and MRC 0156-252, observed 
with HAWK-I in the Y, H and Ks bands. 

MRC 1017-220 (z = 1.768; R.A.: 10:19:49.05, Dec: - 
22:19:58.03, J2000, L 3CHz = 10 28 11 W Hz" 1 ) is our lowest 
redshift targeted HzRG. We note that this HzRG is a broad 
line radio galaxy (Kapahi et al. 1998) and is unresolved both in 
near-infrared and radio dPentericci et alJuOOll) . Investigating ex- 
tremely red objects (EROs; R - K > 6) around high-z AGN, 
ICimatti et al.l (|2000) found an excess of EROs in the close vicin- 
ity of MRC 1017-220 compared to the field. 

We also targeted the field around MRC 0156-252 (z = 2.016; 
R.A.: 01:58:33.63, Dec: -24:59:31. 10,L 3C Hz = 10 2779 W Hz -1 ). 
This HzRG has be en reported to be a qu asar obscured by dust 
dEales & Rawlingslll996ilMcCarthv et al1ll992l) . 



2.2. Observations and data reduction 
2.2.1 . New VLT/HAWK-I data 

The two HzRG fields were imaged between April 2008 and 
August 2009 in Service Mode with VLT/HAWK-I. HAWK-I is 
a wide-field imager on UT4 with a field of view of 7.5' x 7.5' 
equipped with a mosaic of four Hawaii 2RG 2048 x 2048 pixel 
detectors separated by a gap of 15". The pixel scale is 0.1064". 
The field of MRC 1017-220 was observed in Spring 2008 for 
122 min in the Y band (A c = 10210A), 53 min in the H band 
(A c = 16200A) and 33 min in the Ks band (A c = 21460A). 
MRC 0156-252 was also observed in the same set of filters: 
200 min in the Y band (Autumn 2008), 47 min in the H band 
and 33 min in the Ks band (August 2009). In order to avoid 
the gaps between the chips and to have a deeper coverage of 
the immediate surroundings of the HzRGs, MRC 1017-220 and 
MRC 0156-252 were placed near the center of one of the chips. 

Two control fields (hereafter CF1 and CF2, respectively 
centered on R.A.: 11:39:59.66, Dec: -11:24:29.5 and R.A.: 
16:02:06.80, Dec: -17:25:31.7) were also observed from May 
to September 2008 in the Y, H and Ks bands. 
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Table 1. HAWK-I observations 



Field 


R.A." 
J2000 


Dec. a 
J2000 


Band 


Exp.Time 
min 


Seeing 
arc sec 


Limit 
2cr (3cr) 


MRC 1017-220 


10:19:54.18 


-22:18:27.17 


Y 


122 


0.80 


25.91 (25.47) 








H 


53 


0.63 


24.89 (24.45) 








Ks 


33 


0.55 


24.36 (23.92) 


IVIKL- UlJD-ZjZ 


ui.jo.J7. id 


-ZO.UU. J 1 -Ul 


v 


zuu 


n 

U. J^f 


Ifx AO fit QS^ 








H 


47 


0.50 


24.91 (24.47) 








Ks 


33 


0.56 


24.42 (23.98) 


CFl 


11:39:59.66 


-11:24:29.50 


Y 


122 


0.58 


25.80 (25.36) 








H 


53 


0.52 


24.18 (23.74) 


CF2 


16:02:06.80 


-17:25:31.70 


Y 


122 


0.53 


25.78 (25.34) 








H 


53 


0.64 


24.90 (24.46) 








Ks 


33 


0.52 


24.35 (23.91) 



Coordinates are given at the center of the HAWK-I K-band image. 



Part of each F-band observation was observed in pho- 
tometric sky conditions. A standard star, selected from the 
United Kingdom Infrar ed Telescope (UKIRT) faint standards list 
dHawa rden et al.ll2001l) . was observed immediately after the sci- 
ence data to flux calibrate the T-band image. 

Between January and March 2008, the entrance window 
of HAWK-I suffered a degradation and, as a consequence, the 
shadow of the camera spider became visible on the data. The 
data thus contain the convolution of the spider shadow with 
the rotating pupil image, resulting in a cross pattern which re- 
peats the spider symmetry in the backgrouncQ. As advised by 
the HAWK-I User Support Team, we reduced the time for each 
sub-integration to attenuate the cross pattern. However, some 
HAWK-I data taken between April and August 2008 show an in- 
crease in the sky noise independently of the configuration of the 
observations and the cross pattern is therefore very hard to sub- 
tract (e.g., H and Ks-bm&s of CFl). During the data reduction 
phase, we optimized the subtraction of the cross pattern when re- 
moving the background. Unfortunately, the cross pattern in the 
Ks-bwd of CFl could not be properly subtracted and the image 
is unusable. 

The data were r educed using the ESO/MVM (or 'alambic') 
reduction pipelin^ ( Vandame 2004). The image processing fol- 
lowed standard near-infrared reduction steps : dark subtraction, 
division by normalized sky flats, subtraction of the background, 
fringing correction and harmonization of the four chip gains. A 
distortion correc tion was applied t o each chip using stars from 
the USNO-B1.0 dMonet et al.l2003l) catalogue. The images were 
finally stacked using the same astromet ric catalogue (for full de- 
tails on the pipeline, see Vandame 2004]). 

2.2.2. Photometric calibration 

The H and A'i-bands were flux calibrated u sing the 2MASS 
point source catalogue (Skrutski e et alj 120061) for objects with 
11 < K < 14.5 (i.e. 22 stars for MRC 1017-220, 13 for 
MRC 0156-252, 9 for CFl and 37 for CF2). 2MASS stars to- 
tal magnitudes were estimated using SExtractor Kron aperture 
(MAG_AUTO parameter). The derived zeropoints are accurate 
to 0.06 magnitudes. The 2<x and 3<x detection limits of our im- 

1 See http://www.eso.org/observing/dfo/quality/HAWKI/Problems/ 
PupilGhosts.html and http://www.eso.org/sci/facilities/paranal/ 
instruments/hawki/doc/HAWKI-NEWS-2008-07-ll.pdf for details on 
the HAWK-I entrance window problem 

2 http://archive.eso.org/cms/eso-data/data-packages/eso-mvm- 
software-package 




J-K (Vega) 

Fig. 1. Colour-colour diagram (Y - K vs J - K) for the stars from 
a combined UKIDSS/2MASS catalogue (top panel) and for the 
2MASS stars of the MRC 1017-220 field (bottom panel), shown 
as an example for the calibration of the 7-band HAWK-I data. 
Stars are well fitted by an empirical colour-colour relation, (Y - 
K)vega - 1.45 x (7 - K)ve ga (solid line). The standard deviation 
(0.07) of the distribution is shown by the dotted lines. 



ages, determined in randomly distributed 1.5" diameter aper- 
tures, are reported in Table Q] 

We derive the zeropoint of the T-band using the standard 
stars observations and then refine it using empirically derived 
near-infrared colour relations for stars. Using the Wide Field 
Infrared Camera (WFCAM) Science Archivqfl which holds im- 
ages and catalogues of the UKIRT Infrared Deep Sky Surveys 
(UKIDSS), we retrieve stars from a 1 .5x1 .5 square degree region 
(centered around R. A.: 14:00:00, Dec: 10:00:00) and match the 
F-photometry of UKIDSS (data release 5) with 2MASS J and 
^-photometry. We extract stars with accurate photometry in Y, 
J and Ks (less than 0.05 magnitude errors; 1707 stars). Fig. Q] 
(top panel) shows the location of these stars in a (Y - K)y ega vs 



3 http://surveys.roe.ac.uk/wsa/ 
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(J - K) Vega colour-colour diagram. Both UKIDSS and 2MASS 
use Vega photometric systems so Vega magnitudes are used in 
this analysis for consistency. The colour-colour relation for stars 
is well-fit by a simple linear function (Y - K)y eg a - 1.45 X (/ — 
K)vega (standard deviation of 0.07). We extract the 2MASS stars 
from the HzRGs fields and the two control fields and refine the 
zeropoints previously derived from the standard stars observa- 
tions using the above colour-colour relation See the bottom panel 
of Fig. Q] for an example of the calibration method used for the 
MRC 1017-220 field. The offsets applied to the initial zeropoint 
are less than 0.03 mag. We estimate an average 0.06 mag un- 
certainty in the F-band photometry of MRC 1017-220, CF1 and 
CF2 and 0.07 mag for MRC 0156-252. This last field contains 
fewer stars and thus its photometry is slightly more uncertain. 
The 2cr and 3<x limiting magnitudes of the F-band data (deter- 
mined in random 1.5" diameter apertures) are given in TableQ] 

2.2.3. Archival GOODS-S data 

The Southern field of the Great Observat ories Origins Deep 
Survey (GOODS-S; iDickinson et~ai]l2003l) was observed in /, 
H and Ks using VLT/ISAAC from October 1999 to January 
2007. The data we re reduced with the ESO/MVM pipeline by 
the GOODS team dRetzlaffet alJl2010l) . The final data release, 
available since September 2007, includes 24 ISAAC fields in H 
and 26 ISAAC fields in Ks as well as the final H and Ks com- 
bined mosaic^ covering respectively 159.6 and 173.1 arcmin 2 . 

The GOODS- MUlticolor Southern Infrared C atalog 
(GOODS-MUSIC; iGrazian et all l2006bt ISantini et alj [2009) 
is a multiwavelength catalogue of GOODS-S, covering 
143.2 arcmin 2 and cross-correlated optical (u, b, v, i, z from 
Hubble/ ACS and VLT/VIMOS), near-infrared (see above), mid- 
infrared (Spitzer/hac [3.6], [4.5], [5.8], [8.0] and SpitzerfMlPS 
24yum). Spectroscopic redshifts are available for 12% of the 
sources. GOODS-MUSIC provides photometry for 12 of the 
24 fields in H and 22 of the 26 fields in Ks. We make use of 
the GOODS-MUSIC photometry to calibrate the H and Ks 
GOODS-S mosaics. The 2cr and 3cr limiting magnitudes in 
random 1.5" diameter aperture are 25.36 (24.92) and 25.20 
(24.76) for H and Ks. 

As part of the VLT/HAWK-I Science Verification pro- 
grams, two sub-fields of GOODS-S were observed in Y in 
December 2007, centered respectively on R.A.: 03:32:40.92, 
Dec: -27:51:41.6 (355 min; PI: Fontana, A., GOODS 1, here- 
after) and R.A.: 03:32:29.71, Dec: -27:44:38. 6 (145 min; PI: 
Venem ans, B., GOODS2, hereafter). See also Castellan o et alj 
(2010) for further details on the data. The seeing of these im- 
ages is consistent with the one of the ISAAC H and Ks im- 
ages i.e., ~ 0.5 - 0.6". Fig. [2] shows the /fs-band mosaic of 
GOODS-S with the //-band (red) and the two F-band fields, 
GOODS 1 (orange) and GOODS2 (green) overlaid. We reduce 
the data using ESO/MVM pipeline. The images are astromet- 
rically calibrated using a source catalogue extracted from the 
GOODS Ks-bwd mosaic. To flux calibrate the images, we se- 
lect all objects in the field classified as stellar in the NASA/IPAC 
Extragalactic Database (NED). Broad band photometry for these 
stars are taken from the Multiwavelength Surv ey by Yale- 
Chile (MUSYC) catalogues dGawiser et al.ll2006l) . The broad- 
band SED of the stars are fit using stellar templates from the 



4 The reduced single field images and final mosaics are pub- 
licly available at http://archive.eso.org/archive/adp/GOODS/ 
ISAAC_imaging_v2.0/goodsreq.html 

5 Publicly available at http://lbc.mporzio.astro.it/goods/goods.php 
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03:33:00 03:32:48 03:32:36 03:32:24 03:32:12 03:32:00 
R.A. 



Fig. 2. Ks-band combined mosaic of the 26 tiles observed by 
VLT/ISAAC in GOODS-S. The 24 fields also observed in H are 
highlighted by the red dashed line. The two VLT/HAWK-I Y- 
band observed as part of the science verification of the instru- 
ment are overlaid (orange: GOODS 1, green: GOODS2). 

BPGS spectrophotometric atlas dHewett et alj 120061) . For each 
star, the best fitting SED provides the Y magnitude of the star 
and, combined with the flux measured in the image, a zeropoint. 
The image zeropoint is derived from the average of the indi- 
vidual zeropoints and has an uncertainty of 0.05 magnitudes. 
We also independently determine the F-band zeropoint using 
2MASS stars in the GOODS-S field (8 stars) and the relation 
F - K — 1 .45(7 - K) used for photometric calibration in §2.2.2 
and find results consistent within 0.03 mag. The F-band 2<x (3<x) 
limiting magnitudes in random 1.5" diameter aperture are 26.94 
(26.50) and 26.41 (25.97) for GOODS1 and GOODS2 respec- 
tively. 

3. Source extraction 

HAWK-I data, especially those taken prior to May 2009, contain 
crosstalk between the amplifiers of the chips i.e., all the sources 
are repeated on some of the other 32 amplifiers producing a 
series of crater-like artifacts arrayed horizontally dFinger et al.l 
2008). Due to its intensity, only crosstalk produced by the bright- 
est objects is observed out of the background, e.g., in the im- 
age, all sources with K2MASS,vega < 15.5. Regions affected by 
crosstalk were masked before extracting the source catalogues. 
Crosstalk is more prominent in the F-band since the image is the 
deepest and stars are brighter in bluer bands. We therefore first 
identify by eye the crosstalk in the F-band. A map is created 
to flag the crosstalk-affected pixels by 4" x 4" squared masks. 
The flag area accounts for less than 2% of the final F images for 
MRC 1017-220 and CF1 and less than 1% for GOODS-S. The 
CF2 field contains numerous bright stars and ~ 4% of the final 
mosaic is flagged. We also flag regions affected by bright stars 
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that dominate their surroundings as well as the noisy edges of 
the images. 

The source detection was performed using SExtractor 
(Bertin & Arnouts 1996) with a detection threshold of 2<x 
independently for each filter. We used aperture magnitudes 
(SExtractor MAG_APER) within a fixed 2.5" diameter aperture 
to measure colours. Based on the profile of stars in the Ks-band 
data of the HzRGs fields, we estimate that an aperture of 2.5" 
diameter contains about 95% of the source flux. Using such an 
aperture is therefore a good compromise between including as 
much flux from the source as possible, but limiting background 
contamination for faint objects. However, the fraction of a source 
flux contained in a 2.5" aperture strongly depends on the seeing 
of the image and we smoothed the images of different bands to 
the same seeing to ensure accurate colour measurements. 

The 3-band images of MRC 1017-220 and CF2 have sig- 
nificantly different seeings from filter to filter. We therefore 
smoothed our images to the worst seeing. The H and Ks-band 
images of the field of MRC 1017-220 were smoothed to the 0.8" 
seeing of the F-band. Similarly, the F and Ks-band images of 
CF2 were smoothed to the 0.64" seeing of the //-band. We de- 
tected the sources on the unsmoothed images and determined 
the aperture photometry on the smoothed images. We explain in 
§5.3 how we handle upper limits in colours. Total magnitudes 
are determined using SExtractor parameter MAG_AUTO on the 
original unsmoothed images. 

All magnitudes were corrected for Galactic extinction 
(calculated fo r HAW K-I filters) using the dust maps of 
ISchlegel et al] d!998l) and assuming R v = A V /E(B - V) = 
3.1 extinction law of ICardelli et alj (Tl989). All the fields are 
at high galactic latitude (b > 20). Corrections are small for 
the MRC 1017-220 field (0.054, 0.029, 0.019 mag in Y, H 
and Ks respectively), the MRC 0156-252 field (0.013, 0.007, 
0.004 mag) and CF1 (0.033, 0.018, 0.011 mag). Due to a lot 
of dust along the line of sight, corrections for CF2 are, on the 
contrary, rather big (0.344, 0.187, 0.118). We do not apply the 
negligible (< 0.01 in Y and < 0.005 for H and Ks) extinction 
corrections for the GOODS-S field. 

We evaluated the completeness limits of the images using 
a IRAF gallist and mkobjects routines (artdata package) 
to simulate artificial galaxies, both elliptical and spiral galaxies. 
We chose a uniform distribution of galaxy morphologies with a 
minimum galaxy axial ratio b/a of 0.8 and a maximum half flux 
radius of l'.'O. We adopted a de Vaucouleurs and an exponen- 
tial disk surface brightness law for ellipticals and spirals respec- 
tively. For both types of galaxies, we generated catalogues of 
5000 objects and added them to the F, H and Ks images, includ- 
ing Poisson noise. We determined how many artificial sources 
were recovered using the same SExtractor configuration files 
used to detect the real sources. The 90% completeness limits 
for elliptical (spiral) galaxies are 24.0 (23.3), 22.9 (22.2) and 
22.4 (21 .7) in F, H and Ks respectively for MRC 1017-220, 24.4 
(23.8), 23.2 (22.5) and 22.4 (21.8) for MRC 0156-252. 

4. Galaxy number counts in Y, H and Ks 

We derive the differential galaxy number counts in the HzRG 
fields for each band (see Fig. [3). We only consider the deepest 
regions of the final images, i.e., we discard the shallowest re- 
gions resulting from the gap between chips (central 'cross') as 
well as the edges of the image i.e. 12 - 13% of the images. We 
also derive the galaxy number counts combining the two con- 
trol fields: CF1, CF2 and the GOODS-S data for the F and H 
bands. For the Ks band, number counts are determined from CF2 
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Fig. 3. Galaxy number counts in F, H and Ks for the field around 
MRC 1017-220 (red dots), MRC 0156-252 (orange dots) and 
the combined control fields (green diamonds). We do not apply 
completeness correction. We isolate galaxies from stars using 
the stellar index CLASS_STAR in SExtractor. Counts from the 
literature are overplotted for H and Ks. For the F-band, we de- 
rive counts from the two GOODS-S F field and the two control 
fields (CF1 and CF2) until F < 24. Galaxy counts are also de- 
rived from GOODS-S only for sources with 23 < F < 25.5. 
The best fit of the F-band number counts by a two power-law 
model is shown by the dashed line (see §4, also overplotted on 
the //-band number counts for comparison). 
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Table 2. Galaxy number counts in 7-band (N/(mag.deg 2 )) 



Mag. 


1017 


0156 


CFs+GOODS-S 


GOODS-S 


(1) 


(2) 


(3) 


(4) 


(5) 


17.75 


747 


485 


470 


- 


18.25 


896 


1213 


759 


- 


18.75 


1792 


1940 


1228 


- 


19.25 


1792 


3881 


2131 


- 


19.75 


3286 


3759 


3468 


- 


20.25 


5377 


6913 


6141 


- 


20.75 


8812 


7276 


9464 


- 


21.25 


15682 


13461 


13293 




21.75 


18818 


17585 


18531 




22.25 


24494 


26801 


24816 




22.75 


33156 


36746 


36845 




23.25 


31066 


47054 


46742 


44807 


23.75 


35397 


58817 


54652 


54695 


24.25 




78948 




74609 


24.75 








76683 



and GOODS-S, since CF1 was not observed in Ks. We separate 
galaxies from stars using the SExtractor CLASS-STAR param- 
eter. We test the optimal values for this parameter for each im- 
age with stars from 2MASS and USNO in flux bins of 1 mag 
(0.6 <CLASS.STAR< 0.9). 

Table |2]reports (for the first time) the T-band galaxy number 
counts. Counts for the two HzRGs fields (columns 2 and 3) are 
given to the 90% completeness limit. Due to the lower complete- 
ness limits of CF1 and CF2, we first derive counts for Y < 24 
from the four fields (CF1+CF2+GOODS1+GOODS2) i.e., a to- 
tal area of about 200 arcmin 2 (see Table [2] column 4). We also 
derive the counts for 23 < Y < 25.5 from GOODS-S only (col- 
umn 5). 

For the two other bands, w e compare the g al axy counts 
with r e sults from t he lite ra ture: lYan et al.l dl998 | ) : IChen et all 
d2002l): iMov et all d2003h: iMetcalfe et al. i200& for H and 



(2001); Huang et al. 



Maih ara et al.ld2001 | ^urrrmel & Wagnerfd2001l) ; ISaracco et all 



(120011) : lElston et al] d2006l) for Ks. No at- 
tempt is made to correct for the incompleteness, nor correct for 
differences of filter passbands (e.g., K or Ks). 

Number counts have been frequently used to test models 
of galaxy evolution. Although such an analysis is beyond the 
scope of this paper, we note however that the slope shown by 
the F -band counts is consistent with the //-band. Ilmai et all 
( 2007 ) fitted their /-band counts (derived from the AKARI North 
Ecliptic Pole survey) by two power-laws of the form N(mag) = 
aX lQb(ma g -i5) with a ^eak in the slope at J Vega ~ 19.5. Adopting 
a similar fitting function, we find that the T-band counts are well 
fitted by two power-laws: a = 25 + 5; b = 0.45 + 0.02 for Y < 21 
and a = 380 + 5; b = 0.25 + 0.02 for Y > 21. These best- 
fit slopes are plotted in Figf3](top panel) and overplotted on the 
//-band counts (middle panel) for comparison. The Y and H- 
band c ounts are found to be consistent. Similarly to Ilmai et all 
d2007l) . we note that the Ks-bwA number counts show a less 
abrupt change of slope atl9<ZT<21 and cannot easily be 
modeled by two power-laws. 

As shown in Fig. [3] (top panel), the Y number counts of the 
HzRG fields are in good agreement with counts derived from 
the control fields. We note that there is a deficit of sources in 
the MRC 1017-220 field at Y > 23, which cannot be explained 
by a lack of depth of the data since the 90% completeness limit 
is reached at Y = 24 mag. Interestingly, the galaxy counts of 
the targeted HzRG fields show an excess of galaxies in H for 
20.5 < H < 22.5 where we find 27.1 + 0.6% (22.3 + 0.6%) 




1.5 2 
Wavelength (urn) 



Fig. 4. Redshifted (to z = 1.9) spectral energy distribution of 
models of 2 Gyr-old galaxies assuming an exponential declining 
star formation history with r = 0.1, 0.3, 0.5 and 1 (from yellow 
to blue; see Fig. [5] for colours). The transmission curves of the 
HAWK-I filters Y, H and Ks and the position of the restframe 
4000A break (dashed line) are overlaid. 



more sources in the field of MRC 1017-220 (MRC 0156-252) 
than our control fields (CF1 +CF2+GOODS). Similar results are 
found in the Ks-bwd number counts where an excess of sources 
with 19.5 < Ks < 22.5 (27.4 + 0.4%) is found in the field of 
MRC 0156-252 compared to our control fields (CF2+GOODS). 
An excess of sources with 19.5 < Ks < 21.5 (17.8 + 0.6%) 
is also seen in the field of MRC 1017-220. We do not observe 
such excesses in the T-band (which is below the 4000A break at 
Z ~ 2). This suggests that both HzRG fields contain an excess of 
red galaxies. 



5. Search for candidate cluster members at z ~ 2 

5. 1 . YHK colour selection of galaxies atz> 1 .6 

Colour criteria efficiently select high-redshift galaxies in a rela- 
tively narrow redshift range and permit to isolate potential clus- 
ter members associated with HzRGs. At high redshifts (2 < z < 
3), the position of the 4000A break in galaxies falls between the 
J and //-band. Pure near-infrared colours have thus been used 
to isolate candidate protocluster members at this re dshift range, 
e.g. (J — K)ve ga > 2.3, designed by the FIRES team dFranx et al .1 
2003) which allows them to select Distant Red Galaxies (DRGs). 
However, these single colour cuts mostly pick out red passively 
evolving or dusty galaxies and miss galaxies with continuous 
star-formation with little or no dust. 

iKaiisawa et al .1 d2006h defined a two-colour selection tech- 
nique that combined J, H and K colours to select both red, 
passive and blue, star-forming galaxies at z > 2. This criterion 
(/ - K > 2 X (H - K) + 0.5 & J - K > 1.5; Vega system) is 
almost insensitive to dust extinction since the reddening vector 
(E(B - V)) is parallel to the colour selection. 
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Fig. 5. iBruzual & Chariot d2003l) model predictions of different 
stellar populations in a Y—H vs H-K colour-colour diagram (see 
text for details on the models). The dashed line accounts for the 
two-colour selection designed to isolate galaxies at z > 1.6. The 
horizontal line (Y—H > 1.5) shows the single-colour selection to 
separate red passively evolving galaxies (r-YHK galaxies) from 
star-forming ones (b-YHK galaxies). We assume an exponential 
declining star formation history with t = 0.1, 0.3, 0.5, 0.7, 1 
and 5 Gyr (see legend). Top: for increasing redshifts (z — to 
Z = 3, each 0.5 bin marked by coloured points) and a constant 
age of 2 Gyr, the stars corresponding to z — 1.77. The black 
arrow indi cates the reddening v ector E(B — V) = 0.2 as parame- 
terized by ICardelfi et al.l (119891) . Bottom: for various population 
ages (t = 1.5, 2, 2.5, 2.75 Gyr, coloured points with Y—H getting 
redder with ages) at z — 1.77 and z = 2.02. 



At the redshift of the targeted HzRGs, the 4000A break is 
between the Y and the /-band. One can there fore use the Y- 
band i nstead of the /-band in a similar method to Kaiisawa et al. 
(2006) to search for clust er members around MRC 1017-220 and 
MRC 0156-252. We use lBruzual & Charlotl(l2003l) models to de- 
termine the colours of different stellar populations at z — 1 .77 
a nd z = 2.02. In t he models, we assume a solar metallicity and 
a IChabrierl (120031) initial mass function. Fig. |4] shows dust-free 
SED predictions for galaxies with an exponentially declining 
star formation history with r — 0.1 to r — 1 Gyr at z — 1.9. The 
HAWK-I transmission filters (black curves) and the position of 



the 4000A break at z — 1.9 (dashed line) are also overlaid. We 
note that the Y and H filters bracket restframe 4000A. The Y fil- 
ter is narrow (width=0.1/mi, about 3 times narrower than H and 
Ks). The Y - H colour is therefore very sensitive at selecting 
galaxies at the targeted re dshifts. 

Fig.|5]shows the same IBruzual & Chariot) d2003) SED mod- 
els in a Y - H vs H - K colour-colour diagram at various red- 
shifts, with a constant population age of 2 Gyr (top panel), or 
at various population ages (t - 1.5 to 2.75 Gyr) at z — 1 .77 
and z - 2.02 (bottom panel). At z — 0, galaxies have consis- 
tently low Y - H colours (0 < Y - H < 0.3). Beyond z ~ 1.4, 
the Y - H colour becomes redder when the 4000A break en- 
ters the F-band. The H-K colour stays almost constant until 
z ~ 2.5 when the 4000A break enters the H band. By z > 3, the 
Y - H (H - K) colour becomes rapidly bluer (redder) than those 
of galaxies at 1.5 < z < 3. This general trend is observed for all 
models. However, the variations on both the Y - H and H-K 
colours are more pronounced for galaxies with short declining 
star-formation histories i.e. those with stronger 4000A breaks. 
Galaxies with r < 0.5Gyr and z > 1 .6 have Y—H > 1 .5. Galaxies 
with longer star-formation histories never become so red. 

We define a ne w colour criterion, analogous to the 
iKaiisawa et alj d2006f) two-colour JHK selection technique, to 
select galaxies at 1.6 < z < 3 (Fig. [5] dashed line): 



Y - H > 0.7 n Y - H > 2 x (H - K) + 0.3 



(1) 



A single-colour criterion is defined to separate red galaxies 
with old stellar population from blue star-forming galaxies: 



Y-H > 1.5 



(2) 



similar to the single-colour criterion (/ - K)y ega > 2.3 that 
selects DRGs at z > 2. This criterion is shown in Fig. [5] as the 
horizontal solid line. The cut at Y - H = 1.5 was optimized to 
isolate th e reddest popu l ation of galaxies at z > 1.6. By anal- 
ogy with iKodama etail d2007l) notations (r-JHK and b-JHK) 
for galaxies selected by the IKaiisawa et alJ d2006l) near-infrared 
criteria, galaxies selected by equations (l)+(2) are referred to 
as r-YHK galaxies and galaxies selected by equation (1) with 
Y - H < 1.5 are referred to as b-YHK galaxies. 

The YHK criterion is insensitive to dust extinction since it 
has been defined parallel to the reddening vector E(B - V) (see 
Fig.[5J black arrow). However, the dusty star-forming galaxies at 
z > 1 .6 will have similar colours to non dusty, passively evolving 
galaxies, and will also be selected by the r-YHK criterion. Our 
criterion therefore does not enable us to clearly distinguish the 
red passive galaxies from the dusty star-forming ones. 

5.2. Reliability of the YHK colour selection 

We examine the YHK colours of galaxies in GOODS-S us- 
ing the spectroscopic redshifts available in the GOODS-MUSIC 
catalogue (see §2.2.2) and combining the two archival F-band 
(GOODS 1 and GOODS2) images with the H and Ks GOODS- 
S mosaics. We select sources that have a reliable photometry in 
H and Ks (less than 0. 1 magnitude errors), are detected (2<x) in 
the two F-band images, and have a reliable spectroscopic red- 
shift (flag 0: very good or 1: good). Fig. [6] (top panel) shows 
the distribution of these sources in the YHK colour-colour dia- 
gram i.e., 72 stars (stars), 316 galaxies at z < 1.6 (black dots), 
20 galaxies with 1.6 < z < 2.1 (orange circles) and 10 galax- 
ies with z > 2.1 (red circles). AGN (defined as 'BLAGN' or 
'NLAGN' in GOODS-MUSIC) are not considered here. We 
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Fig. 6. YHK colour-colour diagrams of sources with spectro- 
scopic redshifts from the literature (GOODS-S field). The se- 
lection criteria (1) and (2) are shown by the dashed and solid 
horizontal lines respectively. Top: Sources with spectroscopic 
redshifts from the GOODS-MUSIC catalogue (see legend for 
symbols) and stars from the Pickles 1998 st ellar spectra flux li- 
brary. We overplot early-type galaxies from Rett ura et al.l ([2010) 
at z < 1.5 as wel l as th e three passive galaxies at z > 1.8 
from ICimatti et all d2008f) (green and orange triangles respec- 
tively). Bott om: Members of the overdensity at z ~ 1.6 found in 
GOODS-S dKurk et alJ l2009). Colours of the symbols account 
for the spectral class of the group members: red, green and blue 
for early-type, early-type with signs of star formation and late 
type galaxies respectively. 
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Fig. 7. Photometric redshift distribution of the ITiTif-selected 
sources for r-YHK galaxies (top panel) and b-YHK galaxies 
(bottom panel) in GOODS 1+GOODS2 for both the 2cr detec- 
tion limits (solid histogram) and the limit of completeness of the 
HAWK-I data (filled histogram). Photometric redshifts are from 
Santini et al. 2009. 



also overp lot colours of stars from the digital stellar spectra li- 
brary from lPicklesI dl998l) . We note that 13/20 (65%) galaxies at 
1.6 < z < 2.1 are selected by the YHK criteria. The criterion ef- 
ficiently removes foreground galaxies i.e., only 14/316 (< 5%) 
galaxies with z < 1.6 are found in the YHK-selected region. 
4/10 (40%) of galaxies with z > 2.1 are selected by the YHK 
criterion indicating that the selection will be notably contami- 
nated by background objects. 

iRettura et all (1201 0) present a sample of 27 early-type galax- 
ies with 1.09 < z < 1.35 found in the GOODS-S field. 18 
of them are imaged and detected (> 2cr) in the Y band (see 
green triangles in Fig. |6j and have colours consistent with mod- 
els predictions. Only one of them (the highest redshift source at 
z = 1.35) is selected by the YHK criterion showing that con- 
tamination by lower redshift passively evolving galaxies is very 
small. 

ICimatti e~ al. (2008) studied a sample of 13 old, passive 
galaxies at z > 1.4 found in the northern part of GOODS-S 
covered by the Gala xy Mass Assem bly ultra-deep Spectroscopic 
Survey (GMASS; iKurk et al . 1 120081) . Ten of these objects have 
z > 1.6 including three at z > 1.8. The seven others are part of 
an overdensity at z ~ 1.6 presented in Kurk et al. (2009). The 
three passive galaxies at z > 1.8 are well detected in Y (in the 
GOODS2 field), H and Ks (see orange triangles in Fig. [6] top 
panel). Two of them have Y — H > 1.5. The third one is a b-YHK 
galaxy, but has a red Y - H colour, close to the selection limit. 

lKurketai] (T2009) discovered a galaxy overdensity at z — 1 .6, 
with 42 spectroscopically confirmed members in the GMASS 
area. Five galaxies have an early-type spectral class, two galaxies 
are at an intermediate stage (early-type but with sign of star for- 
mation; intermediate type, hereafter) and 35 are late type galax- 
ies. All 42 members are detected in the GOODS2 F-band. Eight 
late-type galaxies have magnitude errors larger than 0.1 mag in 
either Y, H or Ks and will not be considered in our analysis. 
Fig. [6] (bottom panel) shows the colours of the 34 remaining 
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sources in the YHK colour-colour diagram, the colours of the 
symbols indicating their spectral class. 

The five early-type galaxie s and the two ' interm ediate type' 
galaxies are also reported in Cimatt fet al.l d2008l) as passive 
galaxies. All seven sources are selected by the YHK criterion. 
Four early-type and one 'intermediate type' are i-YHK galaxies 
and the last two sources are b-YHK galaxies. 12/27 (44%) of 
the late type galaxies are selected by the criterion showing its 
limitation at z < 1.6. We also note that the [OII]/l3727A doublet 
falls in the F-band at 1.6 < z < 2. If present, and has a large 
equivalent width, the emission line could bias the source F-band 
magnitude and thus its F - H colour, with the galaxy appearing 
bluer than expected. Since the majority (> 80%) of the late-type 
galaxies from Kurk et al. (2009) shows the [Oil] line in their 
spectra, this could explain the blue colours of some of them. 

We further test the YHK criterion using the photometric red- 
shifts (Zphot'i hereafter) available in the GOODS-MUSIC cata- 
logu e. Using the multiwave length photometry of the GOODS-S 
field, iGrazian et al.l d2006bh applied a photometric redshift code 
to their catalogue. They tested their code with the available spec- 
troscopic redshifts a nd found an accuracy of cr z = 0.03 x (1 + z) 
(IGrazian et al.ll2006al) . 

We examine the z P hot of the F/Z/f-selected galaxies within 
the 2<t detection limits as well as within the 90% completeness 
limits of the HAWK-I data i.e., the completeness of passive red 
ellipticals for the i-YHK galaxies and the completeness of spiral 
galaxies for the b-YHK galaxies. Fig. [7] shows the photometric 
redshift distribution of the i-YHK galaxies (top panel) and the 
b-YHK galaxies (bottom panel). We find that 77% (2cr; 71% in 
the completeness limits) of the i-YHK galaxies and 68% (47%) 
of the b-YHK have z p iwt > 1.6 which confirms the efficiency of 
the selection criteria. 

The i-YHK criterion is efficient at selecting galaxies at the 
targeted redshifts with 42% (50%) of the sources with 1.6 < 
Zphot < 2.1. Within the limits of completeness, the b-YHK 
criterion is also efficient with 38% of galaxies having 1.6 < 
Zphot < 2.1. The photometric redshift distribution becomes much 
broader when considering fainter sources and extends towards 
higher redshift sources. At the 2<x limits, 27% of b-YHK sources 
have 1.6 < z < 2.1. We note that if indeed an overdensity of red 
or blue galaxies is present in the surroundings of the HzRGs, the 
sample of F///iT-selected galaxies in these fields would contain a 
higher fraction of sources at the targeted redshifts. The percent- 
ages given earlier are therefore expected to be lower limits. 

5.3. The YHK-selected galaxies 

We apply the YHK criteria to the five fields with YHK cover- 
age i.e., MRC 1017-220, MRC 0156-252, CF2, GOODS 1 and 
GOODS2. We consider sources with a 2<x detection in all three 
bands. The F - H vs H - K colour-colour diagrams for the five 
fields are shown in Fig. [8] Red and blue circles indicate i-YHK 
and b-YHK galaxies respectively. 

We overplot sources detected (within 2cr limits) in H and Ks 
but not in our F-band catalogue. In order to place these sources 
in the colour-colour diagram, we assume a lower limit on the F- 
band magnitude (see arrows, Fig. [8j. We use SExtractor in dual 
mode using the source positions in the Ks-band and deriving 
aperture photometry for these sources on the F-band. For sources 
brighter that our 2cr detection limits but that were not part of 
our F-band catalogue (i.e. beyond our completeness limits), we 
assign the F-band photometry derived from the aperture placed 
at the Ks-band source position. For fainter sources, we assign 



25 



-24:59:30 



35 



■24:59:40 



34.2 01:58:34 33.8 33.6 33.4 33.2 01:58:33 32.8 
R.A. 

Fig. 9. 3-colour image (R, G, B for Ks, H, F) of the 20" X 20" 
field of view around MRC 0156-252 (North is Up, East Left). 
A is the radio galaxy. B and C, previously reported in Pentericci 
et al. 2001, are respectively selected by the i-YHK and b-YHK 
criterion. An additional component, D, is also found aligned with 
MRC 1056-252, B and C but too faint in Ks to be classified by 
our near-infrared criterion. 

them the 2<x limits of our F-band images. Sources with lower 
limits falling in the i-YHK selection area of the colour-colour 
diagram are overplotted in red. These objects are particularly 
interesting since they have very red Y-H colours which suggests 
they have strong 4000 A breaks. 

Within our 2<x limits in F, H and Ks, we find 38 i-YHK 
galaxies (151 b-YHK galaxies) in MRC 1017-220, 105 (191) in 
MRC 0156-252, 47 (196) in GOODS 1, 48 (176) in GOODS2 
and 38 (81) in CF2. We note that since the 2<x magnitude limits 
and area slightly vary from field to field, these numbers are not 
dir ectly compara b le to o ne another. 

ICimatti et all (l2000h looked at the populations of EROs 
((R - K) V ega > 6) in 14 fields around radio-loud AGN at z > 1.5, 
including MRC 1017-220. An excess of EROs was found in the 
field of MRC 1017-220, with three EROs located within 2.5' of 
the HzRG Two of the three were observed in spe ctroscopy in the 
H-band with VLT/ISAAC (ICimatti et all 1 19991) and a 'spectro- 
photometric' redshift {z sp hot) was derived for both sources com- 
bining the ISAAC spectrum continuum and broad-band photom- 
etry. The two sources, J101948-2219.8 (R.A.: 10: 19:47.79,Dec: 
-22:19:46.6, K Vega = 18.7, z sp hat = 1.52 + 0.12) and J101950- 
2220.9 (R.A.: 10:19:49.76, Dec: -22:20:53.9, K Vega = 18.6, 
Zsphot = 1 -50 + 0.25), were both classified as early-type galaxies 
due to their SEDs consistent with no dust extinction. We look at 
the colours of these two objects (see black squares, Fig. [8] bot- 
tom left panel) and both are selected by the b-YHK criterion. 
However, only deeper spectroscopy over a wider wavelength 
range will confirm whether those two targets have redshifts con- 
sist ent with MRC 1 17-22 0. 

iPentericci et all (l200ll) presented near-infrared NICMOS 
imaging of MRC 0156-252. They found two objects within 5" of 
the radio galaxy, the three sources being aligned in the direction 
of the radio axis (i.e. NE/SW). Fig. [9] is a three-colour image 
of the close vicinity of MRC 0156-252. Three objects are ob- 
served near MRC 0156-252. We label the components using the 
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Fig. 8. Colour-colour diagrams Y- H vs H - K for all the fields : Top panels: Control fields: GOODS 1 (left), GOODS2 (middle) and 
CF2 (right). Bottom panels: HzRG fields: MRC 1017-220 (left) and MRC 0156-252 (right). We plot galaxies detected down to the 
2cr detection limits in Y, H and Ks (black dots). i-YHK and b-YHK galaxies are designed by red and blue circles. Arrows indicate 
lower limits for sources detected in H and Ks but not in Y. The radio galaxies are marked by the yellow stars in the bottom panels. 
We also indicate with black open squares the two EROs found in the field of MRC 1017-220 (bottom left panel, Cimatti et al. 1999) 
and the two sources found within 5" of MRC 0156-252 (bottom right panel) corresponding to objects B (the b-YHK galaxy) and C 
(the t-YHK galaxy) in Fig. 



same notation as lPentericci et all d2001l) (A for MRC 0156-252, 
B and C for the two eastern components), adding 'D' for the third 
faint source to the we st of the HzRG. As previously noticed in 
iPentericci et alj d2001l) . component C is much redder than com- 
ponent B. C is a t-YHK galaxy and B is a b-YHK galaxy. B and 
C are indicated by black squares in Fig. [8] bottom right panel. D 
is detected but is very faint in Ks (< 2<x level) and was therefore 
not considered in the candidate selection process. Spectroscopy 
is required to prove the association of these objects with the radio 
galaxy but their physical closeness and colours strongly suggest 
that they are associated with MRC 0156-252. The scale of this 
system is similar t o the structure asso ciated with PKS 1 138-262 
i.e., about 15" (see lMilev et"alll2006l Fig. 2). 



6. Candidates properties 

6. 1 . Surface densities of YHK -selected galaxies 

We compare the surface densities of the TH/T-selected galaxies 
in the five fields to determine whether there is an overdensity 
of sources around the targeted HzRGs. For a direct field to field 
comparison, we cut the catalogues at the completeness limits of 
the shallowest HAWK-I images (see §3; Y < 24, H < 22.9 and 
Ks < 22.4). Table [3] summarizes the densities for both the r- 
YHK and b-YHK galaxies in the five fields assuming Poisson 
errors for the source densities. 

As far as the x-YHK galaxies are concerned, the fields around 
MRC 1017-220, CF2 and GOODS 1 have similar densities. As 
expected, GOODS 2 whi ch contains the overdensity at z = 1.6 
described in iKurk et al.l (120091) is slightly denser than aver- 
age, e.g. by a factor of 1.7 + 0.7 compared to GOODS 1. The 
field around MRC 0156-252 is significantly overdense compared 
to all the other fields: by a factor of 3.1 + 1.1 compared to 
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Table 3. Number densities of YHK-selected galaxies 



Field 


Area 


r-YHK a 


b-YHK* 




(arcmin-) 


(deg 2 ) 


(deg 2 ) 


MRC1017-220 


54.9 


720: 


t220 


2530 


±410 


MRC1017-220 (< IMpc) 


11.0 


1800 


±770 


3240: 


t 1030 


MRCO 156-252 


59.0 


2220 


±370 


2010 


±350 


MRC0156-252 (< IMpc) 


11.0 


2520 


±910 


1800 


±770 


CF2 


54.3 


930 ± 250 


1460 


±310 


GOODS 1 


50.3 


640 ±210 


1570 


±340 


GOODS2 


50.3 


1070 


±280 


1790 


±360 



Within 90% completeness limits; densities were rounded for clarity. 



MRC 1017-220, 3.5 + 1.3 compared to GOODS 1 and 2.1 + 0.6 
compared to GOODS^J 

We also derive densities within 1 Mpc (~ 2', angular separa- 
tion) for both HzRGs corresponding to the classical estimates 
of vir ial radius in the h ighest redshift clusters known to date 
(e.g. iHilton et all 12007). In the close vicinity of MRC 0156- 
252, the i-YHK density is even higher: 3.9 +1.9 times denser 
than GOODS 1 and 2.4 + 1.1 times denser than the (overdense) 
GOODS2 field. 

As far as the b-YHK galaxies are concerned, densities 
are more similar from field to field. However, we find that 
GOODS2 and the HzRGs fields are slightly denser than average. 
MRC 1017-220 is the densest field. The region within 1 Mpc of 
the radio galaxy is 2.2 + 0.8 times denser than CF2. The field of 
MRC 0156-252 is also denser, by a factor of 1.4 ± 0.4 compared 
to CF2. 

We also derive the number counts of the YHK-selected 
galaxies in all the targeted fields. Reference number counts are 
derived for both r-YHK and b-YHK selected galaxies combin- 
ing CF2, GOODS 1, GOODS2 and subtracted from the number 
counts of the targeted HzRGs fields. Fig.fTOlillustrates the num- 
ber counts of this 'excess' of YHK galaxies per 0.5 Ks mag 
bin suspected to be associated with the HzRGs. Errors on both 
the radio galaxies fields and referen ce fields numbe r counts are 
added in quadrature. We assume the iGehrelsl (I 1 9861) small num- 
bers approximation for Poisson distributions. 

As for Table|3] this analysis is made in the limits of com- 
pleteness of the (shallowest) data for a direct field to field com- 
parison. However, due to the selection technique combined with 
our completeness limits, we do not select all YHK galaxies in 
the fields with the selection getting rapidly incomplete at fainter 
magnitudes. To illustrate this, we look at the colours of r-YHK 
galaxies in the field of MRC 0156-252. On average, they have 
(H - Ks) ~ 0.4 and (Y - H) ~ 1.7. The completeness limit 
in Y (our limiting band) is Y — 24 corresponding to H = 22.3 
and therefore Ks = 21.9, despite the data being complete up to 
half a magnitude deeper in Ks. We therefore limit the analysis to 
Ks < 22. 

No significant excess is seen in the i-YHK number counts of 
the MRC 1017-220 field compared to the control fields. The ex- 
cess of b-YHK, already mentioned earlier in this section is also 
observed in Fig. [10] with an excess of Ks > 20.5 sources. The 
field of MRC 0156-252 also shows an excess of blue Ks > 20 
sources. We do not observe an excess of bright red sources 



6 The YHK colour selection is very sensitive to the photometric cali- 
bration of the images and the Y-band zeropoint has relatively large error 
bars (0.07). We note however that even if the zeropoint was offset by - 
0.07, the field would still be overdense in red galaxies by a factor of 2.4 
compared to MRC 1017-220 and 1.6 compared to GOODS2. 
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Fig. 10. Number counts of the excess of YHK-selected galax- 
ies (per 0.5 Ks mag bin) found in the surroundings of 
MRC 1017-220 (top panel) and MRC 0156-252 (bottom 
panel). We first derive reference number counts for the r- 
YHK (red) and b-YHK (blue) galaxies from the control fields: 
CF2+GOODS1+GOODS2 and then subtract them from the 
number counts of the YHK galaxies found in the field of the 
HzRGs. We assume poissonian errors. Densities are given in sq. 
deg. (left axis). Excess number of galaxies for each HzRG field 
(~ 60 sq. arcmin) is shown on the right axis. 



(Ks < 20) in the field of MRC 0156-252 as compared to the 
control fields. The overdensity of red sources in the field of 
MRC 0156-252 becomes prominent at Ks > 21 with an excess 
of 10 to 20 galaxies in the field (by bin of 0.5mag) compared to 
average. At z ~ 2, such magnitudes in Ks f or elliptical galaxies 
correspond to masses of sever al 1O H M dKodama etalj|2007t 
iRocca-Volmerange et alJ20"04l) suggesting that if the overdensity 
detected is indeed associated with MRC 0156-252, the HzRG 
would lie in a structure that already contains very massive, pas- 
sively evolving galaxies. 

6.2. Spatial distribution of the candidates 

The spatial distribution of the if/TT-selected galaxies (detected 
within our 2cr magnitude limits) in the two GOODS fields is 
shown in Fig. [TT] for i-YHK and b-YHK galaxies in red and 
blue large circles. Sources detected in Y, H and Ks (2<x) are 
marked with black dots. Sources detected in H and Ks but not 
in Y, though falling in the i-YHK selection area (i.e. red arrows 
in Fig. [8), are overplotted in small red dots. Due to the deep 
GOODS y-band, there are only two such sources. We observe a 
clear inhomogeneity in the distributions of r-YHK and b-YHK 
galaxies in the northern part of the field (GOODS2). A concen- 
tration of both red and blue galaxies is foun d at the location o f 
the z — 1.6 galaxy overdensity presented in iKurk et al.l (|2009), 
confirming again the efficiency of the YHK criterion at detecting 
z > 1.6 galaxy structures. 

We note that the northern part of GOODS 1 (the south- 
ern part of the GOODS field) is more populated with r-YHK 
galaxies than the rest of GOODS 1. The spatial distribution of 
b-YHK galaxies is more homogenous, but we also observe a 
clear concentration of b-YHK galaxies at the same position as 
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Fig. 11. Spatial distribution of the YHK-selected galaxies in the GOODS fields (GOODS 1+GOODS2). We also indicate with 
op en circles the 6 r- YHK and 12 b-YHK galaxies, spectroscopically confirmed members of the overdensity at z ~ 1.6 reported 
bv lKurk et al.l (l2009h . symbols colours accounting for the object spectral class (see Fig. |6] lower panel). Three b-YHK galaxies 
confirmed at z ~ 1.6 from the ESO GOODS spectroscopy are also indicated by black squares. 



the excess of r-YHK galaxies (i.e. at R.A.~ 53.15 and Dec.~ 
-27 83). We therefo r e confi rm the results of iKurk et all (2009) 
and ICastellano et al.l (120071) that the overdensity spreads south- 
ward. Three of the b-YHK galaxies have a spectroscopically 
confirmed redshift of ~ 1 .6 from the ESO GOODS spectroscopy 
(reported in the GOODS-MUSIC catalogue; see the three black 
squares in Fig. El right panel). 

The spatial distribution of the YHK galaxies around 
MRC 1017-220, MRC 0156-252 and in CF2 is shown in Fig. El 
r-YHK and b-YHK galaxies are shown as red and blue large cir- 
cles (first and second column) respectively. Sources detected in H 
and Ks with lower limits in Y and selected by our r-YHK criteria 
are shown by small red dots. Distances from both HzRGs (yel- 
low star) are indicated on the top and right axis of the four first 
panels. We do not see any specific distribution of YHK-selected 
galaxies in CF2 (Fig. El bottom row). 

For MRC 1017-220, we indicate the two E ROs with spectro- 
photometric redshift from lCimatti et alJ (1 19991) by black squares. 
The r-F/Z/^-selected galaxies have a non uniform spatial distri- 
bution over the field, with a hint of a filamentary distribution in 
the NW-SE direction in which lies the radio galaxy. No clear spa- 
tial inhomogeneity is observed for the b-YHK galaxies except a 
slight excess in the southern part of the field. We note that this 
result is in agreement with previous studies which s how that red 
galax ies are more strong l y clustered than blue ones dDaddi et al.l 
2000: lBrownetlfll200alKong et alJl2006h . 

A visual inspection of the field around MRC 0156-252 con- 
firms the overdensity of red objects around the HzRG, with the 
r-YHK being more concentrated around the HzRG. A concentra- 
tion of b-YHK galaxies is seen in the immediate (r < 2') North- 
West of the HzRG. A second concentration of blue galaxies is 



also observed in the South-East part of the field (at R.A.~29.66- 
29.70, Dec~-25.05). However, no counterpart of this blue galaxy 
excess is seen in the r-YHK galaxy distribution. 

6.3. The colour-magnitude diagram 

The distribution of the candidates in the Y — H vs Ks 
colour-magnitude diagram is shown in Fig. E] Red and blue 
filled circles indicate r-YHK and b-YHK galaxies respectively. 
Candidates in the close vicinity of the HzRGs (r < 2' corre- 
sponding to 1 Mpc at the HzRG redshifts) are shown by larger 
circles in the bottom panels. Both HzRGs are marked by yellow 
starfl We overplot the expected location of a sequence of pas- 
sively evolving galaxies (taken from the Coma cluster at z = 0) 
at z = 1.77 and z — 2.02 for three different formation redshifts 
(zf — 3, 4 and 5). Models at z — 1.6 are also overplotted as the 
GOODS colour-magnitude diagrams (top panels) since a galaxy 
structure at z = 1.6 is known to lie in these fields (see §5.1). 

The right insets of each panel give the histogram of the Y—H 
colours of the candidates galaxies at z > 1.6. It is now well 
known that galaxies in clusters show a strong bimodality in their 
colour distribution, with star-forming non-dusty galaxies being 
on the bluer side, and passively evolving galaxies on the redder 
side (corresponding to the location of the red sequence) sepa- 
rated by a 'green valley' of intermediate type objects. A bimodal 
repartition is suggested in the histograms of the two HzRGs 
(lower panels) with two peaks on both sides of Y — H ~ 1.5 



7 We note that at the redshift of MRC 1017-220 (z = 1.77), the 
[0II]/13727A line is falling in the Y band which can explain the very 
blue Y - H colour of the HzRG. 
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Fig. 12. Spatial distribution of the YH ^-selected galaxies in the fields around MRC 1017-220 and MRC 0156-252 and the control 
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Fig. 13. Colour-magnitude diagrams (Y - H vs K) of the two GOODS fields and the environments of the HzRGs. The r-YHK and 
b-YHK galaxies are marked by red and blue dots. YHK galaxies within 1 Mpc of the HzRGs are highlighted by bigger symbols and 
error bars in the two lower panels. The yellow stars indicate the radio galaxies. The 2cr detection limit of the Ks bands are shown by 
the dotted lines. Models of a red sequence at z — 1 .77 and z = 2.02 are overplotted in the bottom left and right panels respectively 
(solid lines). For information, models of red sequence at z — 1 .6 are also overplotted in the GOODS panels — corresponding to the 
redshift of the galaxy structure confirmed by Kurk et al. (2009) in the GOODS field (dashed lines). Models indicate the predicted 
colours of a passively evolving stellar population (z/ — 3, 4 and 5 from bottom to top). We also show on the right side of each plot, 
the histogram of the Y - H colours of all the candidate cluster members (both red and blue galaxies; open histograms). The filled 
histograms (two lower panels) correspond to T///T-selected galaxies within 1 Mpc of the HzRGs. 



corresponding to the separation criteria that was designed to iso- 
late passively evolving galaxies from star-forming blue ones. 

Half of the red sources within 1 Mpc of MRC 0156-252 have 
colours consistent with red sequence-like galaxies at z — 2.02. 
The scatter of these sources is quite large (over 0.4mag), the in- 
trinsic scatter being enlarged by errors on the F-band photome- 
try, and contamination from background/foreground galaxies. If 
some of the red galaxies re indeed passively evolving galaxies at 
z ~ 2, they would have been formed at very high redshift with a 
formation redshift zj > 3. 



We also note similarities between the colour-magnitude dia- 
gram of MRC 0156- 252 and that of PKS 1 138-262 presented in 
iKodama et aT1(l2007l) . The r-YHK galaxies that lie on the red se- 
quence models in the field of MRC 0156-252 have 21 < K < 23 
with a high fraction of them (~ 40%) lying on the brighter 
end (21 < K < 21.5). This is consistent with the field around 
PKS 1138-262, where the majority of the red near-infrared se- 
lected sources ((J - K)v ega > 2.3) which lay on the red sequence 
had 20.5 < K < 23 and half of them had K < 21 .5. This suggests 
that both fields contain massive red galaxies. 
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At z ~ 2, the 4000A break enters the /-band, /-band 
spectroscopy is thus necessary to prove th at the red galaxies 
are indeed associated with MRC 0156-252. Ikriek et all d2008l) 
present a near-infrared spectroscopic survey of 36 A-bright 
galaxies (Ky ega < 19.7) at z P hot > 2 selected from the MUlti- 
wave l ength Survey by Yale-Chile (MUSYC; iGawiser et"aTl 
2006: lOuadri et al] l2007) and derive spectroscopic redshifts for 
their full sample (1.6 < z spe c < 2.73). Such a spectroscopic cam- 
paign could be conducted similarly in the field of MRC 0156- 
252 since it contains 23 x-YHK (22% of the t-YHK sample) 
and 18 b-YHK galaxies (13% of the b-YHK sample) which have 
A < 21.5. 

iDohertv et al.l (HpiO) describe a spectroscopic campaign in 
the PKS 1 138-262 field. Two DRGs were confirmed to lie at the 
redshift of the radio galaxy (due to the presence of the Ha emis- 
sion line in their spectra). One is a dust-obscured star-forming, 
red, galaxy. The other is an evolved galaxy with little on-going 
star formation. These are the first spectroscopically confirmed 
red galaxies associated with a protocluster at z > 2. These results 
are encouraging for a future near-IR spectroscopic campaign in 
the field of MRC 0156-252. 



7. Summary 

We develop a new purely near-infrared YHK 2-colour selection 
technique to isolate galaxies at z > 1.6 and classify them as 
(i) passively evolving, or dusty star-forming galaxies and (ii) 
star-forming dominated galaxies with little or no dust. We test 
the method using the GOODS-South field, which has been ob- 
served in Y, H and Ks, and for which a large amount of spectro- 
scopic data is available. GOODS-S contains a structure of galax- 
ies at z ~ 1.6 (42 spectroscopically confirmed members so far). 
Applying the near-infrared criteria to the GOODS-S field, we 
recover this structure confirming the efficiency of our new selec- 
tion technique. 

We target the surroundings of two high redshift radio galax- 
ies, MRC 1017-220 (z = 1. 11) and MRC 0156-252 (z = 2.02) 
and a control field using VLT/HAWK-I. The field of MRC 1017- 
220 shows a non-homogeneous filamentary-like spatial arrange- 
ment of red galaxies, and a slight overdensity of blue galaxies. 
MRC 0156-252 lies in an overdense region of both blue and red 
galaxies. This field is 2 - 4 times denser than the other targeted 
fields. 

The red galaxies are clustered around MRC 0156-252 (< 
1 Mpc). The blue galaxies seem to be preferentially distributed in 
two regions; a concentration of b-YHK galaxies is found imme- 
diately at the NW of MRC 0156-252 and another concentration 
is found in the SE part of the field, 4' (~ 2Mpc at z = 2) away 
from the HzRG. Our study of the close vicinity of MRC 0156- 
252 suggests that the radio galaxy has close-by companions with 
three galaxies found aligned with the HzRG (and the radio axis) 
and within 5". This structure is similar to the system found in 
the vicinity of PKS 1 138-262. 

The distribution of the selected red galaxies in a colour- 
magnitude diagram shows that a large fraction of them have 
colours consistent with red sequence-like objects (with zj > 3, 
assuming they contain no dust). The magnitude range of these 
candidates is also similar to the red protocluster members se- 
lected in the vicinity of PKS 1138-262, with red galaxies pref- 
erentially lying on the bright end (K < 21.5). This suggests that 
the systems associated with the HzRGs have already form their 
more massive members. All these results strongly suggest that 
MRC 0156-252 is associated with a galaxy structure at z — 2 



similar to the galaxy system associated with PKS1 138-262 at 
z = 2.16. 
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